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We overview lattice data on d = 1, 2,3 vacuum defects in four-dimensional gluo- 
dynamics. In all the cases defects have total volume which scales in physical units 
(with zero fractal dimension). In case of d = 1, 2 the defects are distinguished by 
ultraviolet divergent non-Abelian action as well. This sensitivity to the ultraviolet 
■ scale allows to derive from the continuum theory strong constraints on the proper- 

f") 1 ties of the defects. The constraints turn to be satisfied by the lattice data. In the 

I 1 5(7(2) case we introduce a classification scheme of the defects which allows to (at 

Oh, least) visualize the defect properties in a simple and unified way. Not-yet-checked 

relation of the defects to the spontaneous chiral symmetry breaking is suggested 
i-C ' by the scheme. 



1. Introduction 

Somehow, it went largely unnoticed that non-perturbative QCD has been 
changing fast. The change is mostly due to results of lattice simulations 
which ask sometimes for novel continuum-theory models, see, in particular, 
. Moreover, in many cases the lattice results arc formulated in specific 
language and do not allow for any immediate interpretation in terms of the 
continuum theory. In this way, there arises a mismatch between richness 
of the lattice data and scarceness of their interpretation in the continuum 
theory. 

An important example of this type are models of confinement. In- 
deed, instantons still dominate thinking on non-perturbative physics on 
the continuum- theory side. On the other hand, it is known from the lattice 
measurements that the instantons do not confine 2 . Moreover, the vac- 
uum fluctuations which are responsible for the confinement have been also 
identified and turn to be monopoles and P-vortices, for review see, e.g., 3 ' 4 . 

There is no understanding whatsoever of these confining fluctuations in 
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terms of the continuum theory. Moreover, if now someone decides to go 
into interpretation of the lattice data there is no regular way to approach 
the problem. The point is that the monopoles and vortices are defined 
on the lattice rather algorithmically, than directly in terms of the gluonic 
fields. The central step is the use of projections which replace, say, the 
original non-Abelian fields by the closest Abelian-field configurations. The 
projection is a highly nonlocal procedure defined only on the lattice and 
blocks out any direct interpretation of the data. 

To circumvent this difficulty, we attempt to summarize here the lattice 
data on the confining vacuum fluctuations entirely in terms of continuum 
theory. Hopefully, this could facilitate appreciation of the results of the lat- 
tice simulations. In particular, we emphasize that the confining fluctuations 
appear to be vacuum defects of dimension lower than d = 4. 

What challenges the continuum theory in the most direct way is a rel- 
atively recent discovery that the monopoles, see 5 and references therein, 
and vortices, see 6 , are associated with ultraviolet divergent non-Abelian 
action. Since gluodynamics is well understood at short distances, this newly 
discovered sensitivity of the vacuum defects to the ultraviolet scale makes 
them subject to strong constraints from the continuum theory 7 . 

The presentation is as follows. In Sect. 2 we summarize lattice data 
on the confining fluctuations. In Sect. 3 constraints from the continuum 
theory are outlined. In Sect. 4 possible relation to dual formulations of the 
Yang-Mills theories is discussed. 

2. Lattice phenomenology 
2.1. Total volume 

Imagine that indeed there exist low-dimension structures in the vacuum 
state of gluodynamics. Which SU(2) invariants could be associated with 
such defects? First of all, we could expect that the total volume of the 
corresponding defects scales in physical units. What this means, is easier 
to explain on particular examples. 

d = defects. In this case, we are discussing density of points in the 
d = 4 space. And the expectation for the total number of the point-like 
defects would be 

Ntot = c A 4 Q CD ■ Vi , (1) 

where V4 is the volume of the lattice and Aqcd can be understood either 
as a position of the pole of the (perturbative) running coupling or, say, 
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as Aqcd = y/&su(2) where <Jsu(2) is the string tension. Appearance of 
Aqcd hi (1) would signal relevance of the fluctuations to the confinement. 

Equation (1) could be readily understood if we were discussing number 
of instantons. However, it might worth mentioning from the very beginning 
that the instantons actually do not belong to the sequence of the vacuum 
fluctuations which we are going to consider. 

d = 1 defects. The d = 1 defects are lines. For the total length one can 
expect 

Ltot = c 1 A 3 QCD -V 4 . (2) 

Such defects can be identified with the percolating monopoles (for lat- 
est data see 8 ' 9 ). Percolation a means that there exists a large cluster of 
monopole trajectories stretching itself through the whole volume of the lat- 
tice. In the limit of infinite volume the percolating cluster also becomes 
infinite. Note also that the monopole trajectories are closed by definition, 
as a reflection of the monopole charge conservation. 

It is worth emphasizing that the scaling law (2) is highly nontrivial and, 
from the point of view of the lattice measurements, represents a spectacu- 
lar phenomenon. Indeed, on the lattice one changes arbitrarily the lattice 
spacing, a while the corresponding bare coupling, g(a) is changed logarith- 
mically, according to the renormgroup. The scaling law (2) implies that the 
probability 9(a) for a given link (actually, on the dual lattice) to belong to 
the percolating cluster is changing as a power of a: 

9(a) ~ (A QCD -a) 3 . (3) 

Thus, in this case powers, and not logs, of the ultraviolet cut off are ob- 
served. In other words, there is no perturbative background to the defects 
which we are discussing and we are addressing directly non-perturbative 
physics. 

d = 2 defects. The defects are now two-dimensional surfaces and for the 
total area the scaling law would read: 

Atot = c 2 A 2 QCD ■ V 4 . (4) 

Such defects can be identified with percolating P-vortices which are 
known to satisfy (4) 4 , for the latest data see 6 . As in all other cases 
discussed here the evidence is pure numerical, though. Because of the 
space considerations we do not discuss here error bars, concentrating only 
on the general picture. Details can be found in the original papers. 

a For theoretical background see, e.g., 10 . 
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d = 3 defects. For a percolating three-dimensional volume we could 
expect: 

= c 3 A QC d- Vi . (5) 
First indications on existence of such defects were obtained recently 11 . 

2.2. Non-Abelian action associated with the defects 

Defects can be distinguished by their non-Abelian action as well. Since we 
have not specified yet the defects dynamically, at first sight, we cannot say 
anything on their action. Surprisingly enough, there exist educated guesses 
concerning the non-Abelian action of the defects based on their dimension 
alone. 

d=l defects. This case is singled out by the consideration that trajec- 
tories correspond to particles. Particles, on the other hand, belong to field 
theory and we may hope to get insight into the properties of the trajectories 
from field theory. And, indeed, the action 

S = M ■ L , (6) 

where L is the length of the trajectory and M is a mass parameter coin- 
cides with the classical action for a free particle of mass M. One may hope, 
therefore, that by evaluating propagation of a particle as a path integral 
with the action (6) one reconstructs the quantum propagator of a free par- 
ticle. And, indeed, this theoretical construction works. Moreover, it is well 
known as the so called polymer approach to field theory, see, in particular, 
12 . Note that the use of the Euclidean (rather than Minkowskian) space 
is actually crucial to evaluate the corresponding path integral. Also, one 
needs to introduce lattice to formulate the theory. 

Although the use of the action (6) does allow to recover the free field 
propagator, the propagating mass turns to be not the same M but is 
equal to 



(const) ( In 7^ 

"prop 



where the constants const, In 7 are of pure geometrical origin and depend on 
the lattice used. In particular, In 7 corresponds to the hypercubic lattice. 

Note that in Eq (7) we reserved now for dependence of the mass pa- 
rameter M(a) on the lattice spacing a. Indeed, tuning of M(a) to ln7/a is 
needed to get an observable mass (7) independent on the lattice spacing a. 
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Thus, our prediction for the action associated with d — 1 defects (which 
are nothing else but the monopole trajectories) is that the action is close 
to 

Smon ~ ' L . (8) 

a 

Indeed, in this way we can explain that the length of the trajectories does 
not depend on a, see Eq (2). 

Prediction (8) does agree with the results of direct measurements of 
the non-Abelian action of the monopolcs 5 . Let us emphasize that the 
prediction (8) does not use anything specific for monopoles and is rooted 
in the standard field theory. Indeed, the polymer approach to field theory 
is no better no worse than other approaches. 

d = 2 defects. There exists simple theoretical argumentation in favor of 
an ultraviolet divergent action of the two-dimensional defects (or vortices) 
as well. Consider the so called gluon condensate 

N 2 - 1 

< (g;,) 2 > « -^4-(! +°M . ( 9 ) 

where G^ v is the non-Abelian field strength tensor and a is the color index. 
Note that the condensate (9) on the lattice is in fact nothing else but the 
average plaquette action 14 . 

The vacuum expectation value (9) diverges as the fourth power of the 
ultraviolet cut off. This divergence is in one-to-one correspondence with the 
divergence in the density of the vacuum energy, well known in the contin- 
uum theory. If one neglects interaction of the gluons, the gluon condensate 
(9) reduces to a sum over energies of the zero-point fluctuations. That is 
why the r.h.s. of (9) in the zero approximation is proportional to the num- 
ber of degrees freedom, that is to the number of gluons. Accounting for the 
gluon interaction brings in perturbative corrections. For details and most 
advanced calculations see the last paper in Ref. 14 . 

Note now that the monopole trajectories with the properties (2) and 
(8) give the following contribution to the the gluon condensate: 

N 2 - 1 

((G£J 2 Un » -^-^(const)(A QCD -a) 3 . (10) 

In other words, monopolcs correspond to a power like correction to the 
perturbative value of the gluon condensate. 

Now, the central point is that there exists rather well developed theory 
of the power corrections, for review see 15 . The leading power-like correction 
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is expected to be associated with the so called ultraviolet renormalon. The 
corresponding contribution is of order 

N 2 - 1 

{{G%) 2 ) uv -ren « -^^(const)(K QCD ■ a) 2 . (11) 

Note that the power of (Aqc d ■ a) here is different from what the monopoles 
give, see Eq (10). 

From this point of view it would be unnatural if the monopoles ex- 
hausted the power corrections. Moreover, if the d = 2 defects with the 
total area satisfying (4) have the action 

Syort = (const) Av J 2 rt , (12) 

then their contribution would fit the ultraviolet renormalon (11). In this 
way one could have predicted (12). The data on the non-Abelian action 
associated with the P-vortices 6 do confirm the estimate (12): 

Svort „ 0.54^ , (13) 

where A vor t is the total area of the vortices. 

d — 3 defects. Proceeding to the three dimensional defects we could 
predict that the story does not repeat itself and there is no ultraviolet 
action associated with the V3. Indeed, if the action would be of order 
Vs/a , the corresponding power correction to the gluon condensate would 
exceed the ultraviolet renormalon, see (11), by (Aqcd ' a ) _1 and contradict 
the theory. 

The data 11 indeed do not indicate any excess of the action associated 
with the V3. One can claim this to be a success of the theory. It is also 
true, however, that the three dimensional defects so far are lacking identity 
in terms of gauge invariant characteristics. 

2.3. Entropy associated with the defects 

The ultraviolet divergences in the action of the monopoles, see (8), and 
of vortices, see (12), suggest, at first sight, that these fluctuations are not 
physical and exist only on scale of the lattice spacing a. Indeed, the prob- 
ability, say, to observe a monopole trajectory of length L is suppressed by 
the action as 

W(L) w exp(— S mon )- (Entropy) ~ exp( — In 7L/a) ■ (Entropy) , (14) 

and in the continuum limit of a — > the suppression due to the action 
is infinitely strong. Since, on the other hand, the observed length of the 
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monopole trajectories does not tend to zero in the limit a — > the suppres- 
sion by action is to be canceled by the same strong enhancement due to the 
entropy. Let us discuss this issue in more detail in cases d — 1,2. 

d=l defects. In this case, the entropy factor is the number N(L) of 
different trajectories of same length L. It is quite obvious that for L fixed 
and a — > 0, N(L) grows exponentially with L 13 . Indeed, trajectory on 
the lattice is a sequence of steps of the length a. The number of steps is 
L/a and at each step one can arbitrarily choose direction. The number of 
directions is determined by the geometry of the lattice. This is the origin 
of the factor ln7/a in Eq (7) valid for the hypercubic lattice. 

Therefore, the observation (8) is nothing else but the statement that 
in case of the monopole trajectories the entropy is fine tuned to the ac- 
tion. Moreover, since there is actually no free parameter in the theory 
(the QCD coupling is running and cannot be tuned) we are dealing rather 
with self tuning of the monopole trajectories. This should be a dynamical 
phenomenon. The author of the present review finds this observation of 
self-tuned objects on the lattice absolutely remarkable. 

d=2 effects. The suppression due to the action is again there, see (13). 
Moreover, from lattice measurements one can conclude that the action is 
not simply proportional to the total area of the vortex but depends also on 
local geometrical properties of the surface, for detail see 16 . 

As for the entropy, the lattice measurements indicate 16 that mutual 
orientation of adjacent plaquettes belonging to the vortices is random. This 
observation, in turn, implies exponential enhancement of the entropy of the 
surfaces. 

Theoretically, one can determine (from indirect argument) entropy of P- 
vortices in case of Z 2 gauge theory In case of the Z 2 theory plaquettes take 
on values ±1 and the P-vortices by definition pierce the negative plaquettes. 
Numerically, 

(Entropy) z 2 -vortices ~ exp( 0.88A tot /a 2 ) , (15) 

where A tot is the total area of the negative plaquettes. 

We see that the action (13) is not tuned to the entropy of ^-vortices 
(15). A straightforward interpretation of this observation is that the en- 
tropy of the P-vortices in SU(2) gluodynamics is lower than that in case 
of the Zi theory. As we will discuss later, the P-vortices in the gluody- 
namics case are populated by magnetic monopolcs. The lower value of 
the entropy could be understood as evidence that at short distances it 
is rather monopoles which are to be considered as fundamental while the 
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P-vortices 'follow' the monopole trajectories. This issue deserves further 
detailed study and discussion. 

2.4. Alignment of geometry and non-Abelian fields 

So far we considered geometry and non-Abelian fields associated with the 
defects separately. The total volume is a geometrical scalar and the total 
action is a scalar constructed on the fields. Of course, observation of invari- 
ant characteristics which mix up geometry and fields would be even more 
interesting. 

d=3 defects. Let us start with d = 3 and assume existence of three 
dimensional volumes to be granted. The volume percolates on the scale 
Aq C£) . However, if we tend a — > then we can think in terms of large 
(in the lattice units) d = 3 defects. What are possible SU(2) invariants 
associated with these volumes? 

In d = 3 the SU{2) gluodynamics is described by three vectors H a 
(a=l,2,3) where the 'magnetic fields' H a are vectors both in the coordinate 
and color spaces. Generically, there are three independent vectors and the 
simplest invariant constructed on these fields is their determinant. However, 
the absolute value of the determinant varies from one point of the d = 3 
volume to another. The invariant which can be associated with the whole 
volume is, obviously, the sign of the determinant b : 

J 3 = sign{e M e abc H?H b k H?} . (16) 

Unfortunately, there are no lattice data which could confirm or reject this 
prediction. 

d = 2 defects. The invariant (16) is uniquely defined as far as all three 
magnetic fields are indeed independent. If there arc only two independent 
vectors the determinant has zero of first order and geometrically we have 
then a closed d — 2 surface, as a boundary of the d = 3 defects. It seems 
natural to speculate that these boundaries are our d — 2 defects, or vortices. 



b The classification scheme of the defects we are proposing below, to our knowledge, has 
not been discussed in literature (an earlier version can be found in Ref. 18 ). However, at 
least partially, it is close to or motivated by well known papers, see, in particular, 17 . It 
is worth emphasizing that we are using the continuum-theory language, assuming fields 
to be continuous functions of the coordinates. Since on the lattice the measurements are 
performed on the scale of the lattice spacing a, the lattice fields fluctuate wildly on the 
same scale. Thus, the underlying assumption is that the continuum-theory language is 
still valid on average, so that the perturbative fluctuations do not interfere with topology. 
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Thus, we come to the prediction that the percolating d = 2 surfaces and 
non-Abelian fields are aligned with each other. In other words, the non- 
Abelian fields associated with the vortices and resulting in the excess of the 
action (12) are predicted to spread over the surface while the perpendicular 
component is to vanish. This prediction works 6 perfectly well. 

d = 1 defects. The next step is an reiteration of the previous one. 
Namely, there could be zeros of the second order of the determinant (16) and 
geometrically zeros of the second order are closed trajectories. Moreover, 
if these closed lines are monopole trajectories (which are indeed closed by 
definition), the non-Abelian field of the monopoles is to be aligned with 
their trajectories. This expectation agrees with the existing data 19 as 
well. 

It is amusing that the monopoles are expected to be locally Abelian. 
Indeed, there is only one independent color magnetic field associated with 
zeros of second order of the determinant (16). Moreover, they are singular, 
see (8). Thus, we are coming to an after-the-fact justification of the use of 
the Abelian projection to detect the monopoles. On the other hand, their 
field is absolutely not spherically symmetrical (which would be the case for 
the Dirac monopoles). And this spatial asymmetry is manifested in the 
measurements 19 . 

2.5. Spontaneous breaking of chiral symmetry 

Note that our classification scheme predicts that the d = 3 defects are 
characterized by invariants which distinguish between left- and right- hand 
coordinates. Of course, on average the regions with determinants (16) pos- 
itive and negative occupy the same volume. (Moreover, in the continuum 
limit the the percolating d — 3 volume, see Eq (5), occupies a vanishing 
part of the whole d = 4 space) . 

As is mentioned above there are no direct measurements of the sign of 
the determinant (16) within the d=3 defects. However, it is known that 
removal of the P- vortices restores chiral symmetry 20 . Also, there exists 
independent evidence that some 3d volumes are related to the spontaneous 
symmetry breaking in case of SU(3) gluodynamics 21 c . However, it seems 
too early to speculate that the two 3d volumes in question are actually the 
same. 



c An obvious generalization of our invariant (16) to the SU(3) case would be 
sign {fab^ ikl HfH b k Hf}. 
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3. Theoretical constraints 

3.1. Consistency with the asymptotic freedom 

It is worth emphasizing that there is no developed theory of the defects 
considered. So far we summarized lattice observations. Moreover, we have 
been even avoiding discussion how the defects are defined and observed on 
the lattice, for reviews see, e.g., 3 ' 4 * 22 . Let us only mention that the guiding 
principle to define the defects was the search for effective infrared degrees 
of freedom responsible for the confinement. 

However, what is most amusing from the theoretical point of view is 
that the defects have highly non-trivial properties in the ultraviolet. The 
ultraviolet divergence in the action, see (8), (13) are most remarkable. 

Indeed, Yang-Mills theories are well understood at short distances. The 
only divergence which is allowed on the fundamental level is that one in 
the coupling a s . Thus, one could argue that all the ultraviolet divergences 
are calculable perturbatively in asymptotically free theories. Moreover, 
this statement seems rather trivial. What is actually not so trivial is that 
on the lattice one can also consider power-like divergences, see, e.g., (10), 
(11). In the continuum theory, power- like ultraviolet divergences usually 
are used, at best, for estimates. On the lattice, the ultraviolet cut off is 
introduced explicitly and one can treat power-like divergent observables in 
a fully quantitative way, see, e.g., (9). This extends in fact the predictive 
power of the theory. 

It is, therefore, no surprise that using the asymptotic freedom one can 
derive strong constraints on the properties of the vortices 7 . 

3.2. Classical condensate < 4> magn > 

Let us consider first the vacuum condensate of the magnetically charged 
field < 0\4>magn\0 >• Of course, in the YM theory there is no fundamental 
magnetically charged field. However, the monopole trajectories are ob- 
served on the lattice. Using the polymer approach to field theory we can 
translate the lattice data on trajectories into a field-theoretic language, see 
Ref 23,7 and references therein. The only assumption is that there exists 
an effective field theory for the magnetically charge field. 

In particular, the percolating cluster corresponds to the classical vacuum 
expectation value <(/)>. One can derive 23 : 

<4>magn> 2 ~ \p V erc ~ (const) Aq CD (a ■ Aq C d) ■ (17) 
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3.3. Vacuum expectation value < \4> m agn\ 2 > 

One can also expect that there exist quantum fluctuations. And indeed, 
apart from the percolating cluster, there observed finite monopole clusters. 
Which arc naturally identified with the quantum fluctuations. A basic 
characteristic for these clusters is again their total length. By definition: 

L to t = L perc + Lfi n = Pp er c • V\ + Pfin ' V4 , (18) 

where p per c,Pfin are called the densities of the percolating and finite 
monopole clusters, respectively. 

Using the polymer approach to the field theory one can express the 
vacuum expectation value of the magnetically charged field in terms of the 
monopole trajectories 23 : 

< \<f>\ 2 > = ^(Pperc + Pfin) ■ (19) 

Instead of deriving this relation (which is also quite straightforward) let us 
explain why (19) is natural. Concentrate on the quantum fluctuations, that 
is on pfi n . Moreover, consider small clusters L « A.q Cd . Then there is no 
mass parameter at all and on pure dimensional ground one would expect: 
/ const const 

\Pfin) dimension = a d-l = fl 3 ( 2 ™ 

where d is the number of dimensions of space and we consider the d = 4 
case. If this dimensional estimate held, then the vacuum expectation value 
would be quadratically divergent in the ultraviolet: 

< \4 > \ ^dimension ~ a ■ (21) 

And we rederive the standard quadratic divergence in the vacuum expec- 
tation value of a scalar field. 

Now, the central point is that although we call all these estimates 'di- 
mensional' or 'natural' we are not allowed to have (21). Indeed, Eq. (21) 
would hold for an elementary scalar field. However, we are not allowed to 
have new elementary particles at short distances. Because of the asymp- 
totic freedom, there are only free gluons at short distances. And what we 
are allowed to have for the vacuum expectation value in point? Clearly: 

< |^| 2 >allowed ~ ^QCD ■ ( 22 ) 

In terms of the monopole trajectories (which are our observables) Eq (22) 
reduces to: 

Pfin ~ • (23) 
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It is most remarkable that the data 8,9 do comply with (23)! 

3.4. Branes 

It is of course very gratifying that the data comply with the constraint 
(22). On the other hand, the reader may feel that our summary of the 
phenomenology looks self contradictory. Indeed, first we observed that the 
non-Abelian monopole action corresponds to a point-like particle, see (8). 
But then we said that there should be no new particles, and the data agree 
with that constraint. 

Still, there is no contradiction between these observations. Rather, 
taken together they amount to observation of a new object, which can 
be called branes. Indeed, cancellation of In 7 /a in the equation for the mass 
is needed to balance the entropy at very short distances of order a. Eq (22), 
on the other hand, is a global constraint. The geometrical meaning of this 
constraint is actually transparent. Namely, it means that on large scale the 
monopoles live not on the whole d = 4 space but on its d = 2 subspace. 

Within the classification scheme of the defects which we discussed above, 
this association of the monopoles with surfaces is automatic. Algorithmi- 
cally, however, the trajectories and surfaces are defined independently And 
the fact that the monopole trajectories do belong to surfaces is highly non- 
trivial from the observational point of view. 

Thus, what is observed on the lattice are d = 2 surfaces populated 
with "particles" (better to say with the tachyonic mode of the monopole 
field). When we call this objects branes we emphasize that the affinity of 
the monopoles to the surfaces remains true even at the scale a. Traditional 
discussions of the P-vortices and monopoles emphasize, on the other hand, 
physics in infrared, and one talks about 'thick vortices', see, e.g. 4 . 

3.5. Implications for models 

Abelian Higgs model. The lattice data on monopoles are usually interpreted 
in terms of an effective Abelian Higgs model, see in particular 24 and ref- 
erences therein. Our Eq (17) implies, however, 

< <f>magn > ~ (a ■ KqCdY' 2 ^QCD (24) 

It is most remarkable that the classical condensate vanishes in the contin- 
uum limit a — > 0. Nevertheless, the heavy quark potential at large distances 
generated by the monopoles remains the same since it is determined entirely 
by Pperc which scales in the physical units! 
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To establish a relation to the standard fit to the Abelian Higgs model 
one should use, most probably, matching of the two approaches at some a 

d 

Gauge invariant condensate of dimension two. Eqs (22), (23) provide 
us with a value for a gauge invariant condensate of dimension two: 

< \4>magn\ 2 > ~ d-Pfin ~ ■ ( 25 ) 

In terms of the fundamental variables, condensate of dimension two was 
introduced 25 as the minimum value along the gauge orbit of the gauge 
potential squared (^)m in ■ While the < {A^) 2 > m i n is contaminated 
with perturbative divergences, the condensate (25) provides us, on the phe- 
nomenological level, directly with a non-perturbative condensate of dimen- 
sion two. Existence of a gauge invariant condensate of dimension two is 
crucial for models of hadrons, see, in particular, 26 . On theoretical side, 
the nonperturbative part of the condensate of dimension two in the Hamil- 
tonian picture is related to the Gribov horizon 27 . 

Percolation. The fact that the monopole action is ultraviolet divergent, 
see (8), allows to consider them as point-like and apply percolation the- 
ory. In particular, one derives in this way the spectrum of finite monopole 
clusters as function of their length Lf in : 

N(L fm ) ~ j±- , (26) 

fin 

which agrees with the data, for details see 23 . 

Moreover, it is a common feature of percolating systems that near the 
phase transition, in the supercritical phase, the density of the percolating 
cluster is small. Then the observation (3) can be interpreted as indication 
that in the limit a^Owe hit the point of second-order phase transition. 

Note also that second order phase transition is associated usually with 
a massless excitation. At first sight, there is no massless excitation in 
our case, however. The resolution of the paradox is that percolation (or 
randomization) happens now on the ultraviolet scale. Respectively, all the 
masses are to be measured in the lattice units, 1/a. In this sense a glueball 
mass, 

mllueball ~ ^2 " ( fl " A QCD) 2 , (27) 



The author would like to acknowledge discussion of this point with M.N. Chcrnodub 
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in the limit a — > corresponds to a massless excitation. Moreover, one 
expects that the glueball mass controls the scale of fluctuations of the 
monopole clusters. But these expectations have not been checked yet (see, 
however, 2S ). 

Singular fields. Although we are repeatedly emphasizing that the branes 
are associated with singular fields it should also be mentioned that the non- 
Abelian fields considered are significantly smaller than the corresponding 
projected fields. In particular, the Dirac monopole would have a larger 
action: 

S Dir ~ \ ■ - > In 7- (28) 
g z a a 

where g 2 is the gauge coupling and we consider the limit g 2 — ► 0. Already 
at presently available g 2 the action calculated in terms of the projected 
Abelian fields (corresponding to the Dirac monopole) is a few times larger 
than the actual non-Abelian action which determines the dynamics and 
which we discussed so far. The same is true for the vortices. 

This distinction - in terms of the action - between the Dirac and lattice 
monopoles is very important from theoretical point of view. The observed 
monopoles are associated with singular non-Abelian fields but these fields 
are no more singular than ordinary zero-point fluctuations, or perturbative 
fields. The Dirac monopoles, on the other hand, in the limit of g 2 — > would 
be more singular than the perturbative fields. According to the standard 
ideas of the lattice theories, such fields could actually be removed without 
affecting the basic physical content of the theory. The lattice monopoles 
and vortices, on the other hand, cannot be removed without affecting per- 
turbative fluctuations as well. In the next two sections we will consider this 
issue in more detail. 

4. Towards duality 

It is a well known that topological excitations of a 'direct' formulation of 
a theory may become fundamental variables of a dual formulation of the 
same theory. Examples can be found, e.g., in the review 29 . Little, if 
anything, known theoretically on the dual formulation of the Yang-Mills 
theories without supersymmetry. Nevertheless, generically one might think 
in terms of branes 30 . In case of supersymmetric extensions of YM theories 
the branes are classical solutions. One could speculate that the branes 
discussed in this review are 'quantum branes'. Of course, it remains a pure 
speculation until something definite could be said on the properties of the 
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quantum branes on the theoretical side. It is amusing, however, that there 
is a sign of duality between the branes discussed in the preceding section 
and high orders of perturbation theory 31 . 

5. Long perturbative series 
5.1. Expectations 

Let us start with reminding the reader some basic facts about perturba- 
tive expansions, for detailed reviews see, e.g., 15 . A generic perturbative 
expansion for a matrix element of a local operator looks as: 

oo 

( O ) = (parton model) ■ (l + ^ a n a n 8 ) , (29) 

71 = 1 

where we normalized the anomalous dimension of the operator O to zero 
and a s is small, q s <1. 

In fact, expansions (29) are only formal since the coefficients a n grow 
factorially at large n: 

Kl ~ c?-nl , (30) 

where Cj are constants. Moreover, there are a few sources of the growth 
(30) and, respectively, Cj can take on various values. The factorial growth of 
a n implies that the expansion (29) is asymptotic at best. Which means, in 
turn, that (29) cannot approximate a physical quantity to accuracy better 
than 

/ \bo/ci 

A ~ exp(-l/ Ci a s ) ~ [A 2 QCD -a 2 j , (31) 

where 60 is the first coefficient in the /^-function. To compensate for these 
intrinsic uncertainties one modifies the original expansion (29) by adding 
the corresponding power corrections with unknown coefficients. 

In case of the gluon condensate the theoretical expectations can be 
summarized as: 



(0| -J-^_( G ^) 2 |0} « a / \ ' (1+ ]T a n a n s + (const)a 4 ■ A 4 QC D ) , 

n— 1 

(32) 



where 



b a s 



and terms proportional to Aq C£) correspond to < 0|(G^) 2 |0 > so ft which 
enters the QCD sum rules, for review see 32 . 
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A conspicuous feature of the prediction (32) is the absence of a quadratic 
correction, compare (11). Thus, we are seemingly coming to a contradiction 
between the lattice branes and continuum-theory perturbation theory. Let 
us, however, have a closer look at the problem. 

5.2. Numerical results 

Numerically, this perturbative expansion for the gluon condensate was stud- 
ied in greatest detail in the papers in Ref 14 (especially in the latest one). 
The results can be summarized in the following way. On the lattice, the 
gluon condensate is nothing else but the average plaquette action. Repre- 
sent the plaquette action (P) as: 

a A {P) « P p N ert + b N a 2 A 2 QCD + c N a A k% CD , (33) 

where the average plaquette action (P) is measurable directly on the lat- 
tice and is known to high accuracy, Pp ert is the perturbative contribution 
calculated up to order N: 

n—N 

p p N ert = i - Yl p-9 2n > ( 34 ) 

n=l 

and, finally coefficients &at,c/v are fitting parameters whose value depends 
on the number of loops N. Moreover, the form of the fitting function 
(33) is rather suggested by the data than imposed because of theoretical 
considerations. 

The conclusion is that up to ten loops, N = 10 it is the quadratic 
correction which is seen on the plots while cm are consistent with zero. 
However, the value of 6jy decreases monotonically with growing N 14 . The 
factorial divergence (30) is not seen yet and perturbative series reproduces 
the measured plaquette action at the level of 10~ 3 . Finally, at the level 
10~ 4 the h-QCD term seems to emerge 14 . 

5.3. Implications 

Thus, there is a fundamental difference between the instantons and branes. 
The instantons correspond to the 'soft' gluon condensate and are hidden 
in the (Aqcd • a) 4 corrections which are not calculable perturbatively. In 
short, instantons are added to perturbation theory. 

The branes, on the other hand, appear to be dual to long perturbative 
series. If one is able to calculate many orders of perturbation theory, there 
is no need to account for the branes as far as local quantities are concerned. 
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This might first sound disappointing for those who is beginning to be- 
lieve in the important role of the lattice branes. In fact it is not disappoint- 
ing at all. To the contrary, we have first firm piece of evidence that the 
branes belong to a dual world. To reiterate: the instantons belong to the 
'direct' formulation and they are added to the perturbation theory of the 
direct formulation. Branes belong to the dual formulation. Adding them 
to the perturbation theory of the direct formulation would be mixing up 
two different (dual to each other) representations of the same theory. 

Actually, the very existence of fine tuned branes could not be under- 
stood within the direct formulation of the YM theory. However, if there 
were an existence theorem for a dual formulation, the fine (or self-) tuning 
was implied by the theorem. Reversing the argument, we can say that ob- 
servation of the fine tuning might indicate existence of a dual formulation. 
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